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rr^ Abstract: We consider a model for the Higgs sector with two scalar doublets and a 

,^ softly broken Z2 symmetry, the Stealth Doublet Model. This model can be seen as a 

generalization of the Inert Doublet Model. One of the doublets is the Higgs doublet that 

'T^ participates in electroweak symmetry breaking and couples to fermions. The other doublet 

Tjj- does not couple to fermions at tree level and does not acquire a vacuum expectation value. 



The broken Z2 symmetry leads to interesting phenomenology such as mixing between the 

two doublets and charged and CP-odd scalars that can be light and have unusual decay 
T^ channels. We present theoretical and experimental constraints on the model and consider 

^^ the recent observation of a Higgs boson at the LHC. The data on the -?/ — )• 77 channel can 

be naturally accommodated in the model, with either the lightest or the heaviest CP-even 

scalar playing the role of the observed particle. 



1 Introduction 

The ATLAS [1] and CMS [2] experiments at the Large Hadron Colhder (LHC) have dis- 
covered a new particle that exhibits all the features of a Higgs boson (for the most recent 
data see e.g. [3-5]). It will require hard work to uncover if this is the Standard Model 
(SM) Higgs boson or not, but the data seem to be compatible with the SM. If it is a Higgs 
boson, but not the one predicted by the SM, then there is likely an extended Higgs sector 
with additional Higgs bosons. 

Arguably, the most "standard" Higgs scenarios beyond the SM are the minimal su- 
persymmetric Standard Model (MSSM), general two-Higgs doublet models (2HDM), or 
perhaps the next-to- minimal supersymmetric standard model (NMSSM); see in particular 
ref. [6] for a recent review of 2HDMs. These models all predict a similar set of addi- 
tional Higgs bosons, with associated dominating production and decay channels, and most 
searches are devoted to these channels. 

There is a real possibility that Nature is not described by one of these standard scenar- 
ios, and if it is not, the standard searches may not be appropriate. It is therefore important 
to consider alternative scenarios. One such alternative is the Inert Doublet Model (IDM) 
[7-9], where there is a conserved Z2 parity, such that while one scalar plays the role of a 
SM-like Higgs boson, the other scalars do not couple to fermions. In particular, the lightest 
scalar is a stable dark matter candidate. 

In this paper we present the Stealth Doublet Model (SDM), a generalization of the 
IDM where the TLi symmetry is softly broken, and discuss its impact on Higgs physics at 
the LHC. In particular we study the SDM in relation to the LHC discovery and exclusion 
results. When the Z2 symmetry of the IDM is broken, the lightest scalar is not stable, and 
couplings of the fermiophobic particles to fermions are generated at one-loop level. The 
resulting phenomenology is very different from the standard scenarios discussed above, and 
as we will show below is also compatible with the discovered particle at LHC. 

The Z2 symmetry is broken softly by a dimension-two term in the scalar potential. We 
consider this as a way to parametrize our ignorance of the symmetry breaking mechanism, 
in a similar way to how soft supersymmetry breaking is parametrized in e.g. the MSSM. 

The particle content of the SDM is the same as in CP-conserving 2HDMs: there are 
two CP-even neutral scalars, l^ and H^ ^ one of which should be the Higgs boson discovered 
at the LHC, a CP-odd neutral scalar A'^, and a charged scalar H . The interactions of 
these particles are quite different from those of 2HDMs, however. 

In our model, the J^ and K have no tree-level couplings to fermions: they are fermio- 
phobic. All their couplings to fermions are instead generated at the one-loop level. Their 
dominating production and decay modes can therefore be different than in the standard 
scenarios. Consequently, A^ and H can be lighter in the SDM than in standard scenarios, 
since flavor constraints and LEP limits do not apply. For example, if the charged scalar 
is the lightest scalar, its main decay is typically H^ — )■ W^j. In addition, electroweak 
precision tests (EWPT) allow both lighter or heavier hP,H^ in our model than in the SM. 

In the rest of this paper we describe the SDM and its parameters, and we outline 
the requirements for soft symmetry breaking. We consider the various constraints on the 



model and the observed Higgs boson signal at LHC. Finally we briefly consider the charged 
scalar and its decay and production channels. We leave detailed analyses of the model for 
upcoming papers [10], where we will also compute all decay widths of the fermiophobic 
scalars and discuss LHC phenomenology in more detail. A preliminary presentation of this 
model can be found in [11]. 

2 The Stealth Doublet Model 

In brief, the model consists of adding another doublet to the Standard Model scalar sector, 
and introducing a discrete Z2 parity between the two doublets. This parity is broken 
softly, which imposes certain relations on the parameters of the scalar potential. The scalar 
potential is thus the potential of two- Higgs doublet models (2HDM) with two hypercharge 
Y = 1 scalar doublet fields $i = {(l)f , cp'-)'^ , 

V = mfi$|$i + m22<l>2*^2 - ["ii2*^i^2 + h.c] 

+ ^Ai($l$i)2 + '^\2i^l^2? + A3($I^l)(^5$2) + A4($I$2)($l$l) 

+ UA5(^I$2)'+[A6(^I$i) + A7($5^2)]^I^2 + h.c.| , (2.1) 

where we only consider CP-conserving models and therefore take all parameters to be real. 
The Z2 symmetry transformation can be taken as 'I'l — t- <1>i and <1>2 -^ —^2- This parity is 
broken by the last three "hard-breaking" terms in eq. (2.1) and by the soft-breaking term 
ml2^i^2 + h.c. 

Furthermore, in this model the second doublet does not get a vacuum expectation 
value (vev) and is therefore not really a Higgs doublet. However, because the Z2 symmetry 
is broken the two doublets can mix. Studying a model where the vev resides solely in one 
of the doublets is equivalent to working in the Higgs basis (see e.g. [6, 12, 13]). The Higgs 
basis is often useful for analyzing a general two-Higgs doublet model, but in our model, 
this is the physical basis, with v = vi ^ 246 GeV. Similarly to the IDM, the SDM does 
not have a parameter tan/3 = v^jvx^ and cannot be obtained by simply taking the limit 
tan/3 — )• or tan/3 — )• 00. 

Minimizing the potential yields the conditions 

and ?n,22 is thus a free parameter. Further constraints on the parameters will be discussed 
below. 

Generally in two-Higgs doublet models where CP is conserved, there are two CP-even 
neutral scalars, one CP-odd neutral scalar J^ ^ and a charged scalar H . The two doublets 
can be written in unitary gauge as 



^1 = U + 0? ' ^2 = U^ + i^° , (2.3) 





where (/>5 2 are the neutral CP-even interaction eigenstates, whose mass matrix is not diag- 
onal: 

M^-i'"': /-; yj'^i y- x (2.4, 

where A345 = A3 + A4 + A5. The physical CP-even mass eigenstates are given by 

H^ = (j)\cosa + <i)ls\na (2.5) 

h^ = -(t)\sma + (tP2Cosa, (2.6) 

where the mixing angle a is obtained by diagonalizing the matrix A^^. The couplings of 
h^ and H^ to Z^Z^ and W W^ are thus suppressed compared to the SM Higgs by factors 
sin a and cos a, respectively. Taking H^ to be the heavier state, the masses of h^, H^ are 

ml = cljn\ + sl^v^Xi + c^w^As - 2saCav'^\Q (2.7) 

m\j = s1jn\ + c\v^\i + s^f ^As -|- 2saCaV^\, (2.8) 

where we defined the abbreviations Sa = sino, Ca = cos a. The masses of the remaining 
states are 

^A = mjj± - ^v\X5 - A4) (2.9) 

"^H± = "i22 + 2^^^3. (2.10) 

One may solve eqs. (2.7-2.10) for the parameters Ai^3^4^5 so that the masses of the scalars 
can be used as model parameters. The mixing angle is given by 

sin 2a = f'^\ . (2.11) 

mj^-ml 

From eqs. (2.4) and (2.11), we see that when the Z2 parity is exact (i.e. Xq = 0, to^2 = 0), 
there is no mixing of the CP-even states and we recover the Inert Doublet Model. In fact, 
all our results reduce to the IDM when letting Ag — )■ and either sin a — )■ 1 , cos a — )■ or 
sin a —7- 0, cos a — >• 1. In this sense, our model is a generalization of the IDM. 

The Z2 symmetry, if exact, would forbid flavor- changing neutral currents (FCNC). If 
the symmetry is broken, there could potentially be large FCNC. In this model the symmetry 
is only broken softly, i.e., the symmetry is restored at very high energies compared to the 
soft-breaking mass scale |mi2|, and FCNC only occur at two- loop level [10]. 

The requirement of soft breaking introduces additional constraints on the parameters, 
eq. (2.2) shows that the soft-breaking mass scale is proportional to the hard breaking 
coupling Ag, so one might naively expect that if there is soft breaking, there is also hard 
breaking. This is not correct, however, due to the freedom to choose a basis for the 
doublets, i.e., to choose vi and V2 [12-14]. There may exist a basis in which the hard 
breaking parameters Ag = A7 = while 771^2 7^ such that the symmetry is broken softly. 
If this is the case, then the symmetry is broken softly in any basis. 



We analyze this situation by using the methods of refs. [12-14] and find that the 
conditions for soft Z2 breaking, which hold if Ag 7^ 0, are 

(Ai - A2) [A345(A6 + A7) - A2A6 - A1A7] - 2(A6 - A7)(A6 + A7)' = 0, (2.12) 

(Ai - \2)ml^ + (As + h){ml^ - ml^) / 0. (2.13) 

These equations imply that since the parameters are real, there exists a maximum value 
^max £qj, ^^ ^j^ couplings Aj which fulfill eqs. (2.12, 2.13) break the Z2 symmetry softly. 
These relations are thus part of the specification of the model, and in the following we 
always choose parameters to respect this requirement. 

The scalar potential (2.1) has ten free parameters when requiring all couplings to be 
real. Two are eliminated by the minimization conditions (2.2) and one by the conditions 
for soft breaking (2.12, 2.13). We will choose the masses of the scalar states rrih, rrin, itla 
and mu± as four of the remaining parameters. To specify the amount of Z2 breaking, we 
choose to use sin a, which is related to Ae through eq. (2.11). Two more parameters need 
to be specified. We take these to be A3 and A7, since A3 determines the coupling between 
h^/H^ and H^. In summary, we will use as the seven parameters of the model 

mh, niH, ruA, rnH±, sin a, A3, A7. (2-14) 

The simplest solution to (2.12) is to choose A7 = Xq and A2 = Ai. The inequality 
(2.13) then becomes 771,22 ¥^ "^11 • ^^ ^^^ numerical calculations we will often choose A3 
to take the values A3 = 0, 2m'j^±/v'^ and Am'jj±/v'^, corresponding to m^g = fnj^^, and 
— m^±, respectively. We will also vary A3 and A7, within theoretically allowed regions, to 
deduce their impact on the signal strengths for h^ /H^ — )• 77. In ref. [10], we will study the 
solutions of eqs. (2.12, 2.13) in more detail. 

As a final step we must specify the Yukawa couplings of the model. The most general 
Yukawa Lagrangian in the Higgs Basis reads, in terms of the fermion flavor eigenstate [13] 

-£Yuk = KqLl^iEr + k^Ql^i Ur + K^Ql^iDr 

+ P^Ll^2Er + p'iQL^2 Ur + P^Ql^2Dr (2.15) 

where $i = — io"2^i • -'■'^ order to obtain mass eigenstates, the matrix Kq can be diagonalized 
by a biunitary transformation to obtain the diagonal mass matrix M^ for fermions F = 
U, D, L. The correspondingly transformed p^ matrix will in general be non-diagonal and 
will generate FCNC. To avoid FCNC at tree level, we impose positive Z2 parities on the 
fermions. This enforces p = at tree level. The doublet $2 is therefore fermiophobic, and 
the states H , A^ , and (f)^ do not interact with fermions at tree level. Fermions acquire 
mass through Yukawa couplings with the Higgs doublet <I>i only. In this sense, our model 
is a Type I 2IIDM. The Higgs Yukawa Lagrangian is then 

-CYnk = ^-i^fi'fiCaH - Sah) (2.16) 

for all fermions /. The soft breaking terms will lead to couplings between <I'2 and fermions 
at one-loop level. The resulting p^ matrices are diagonal at one- loop level [10]. At higher 
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Figure 1. Examples of allowed regions in parameter space taking into account theoretical con- 
straints and S and T values as a function of to/j± and ttia- The color displays the deviation from 
the center of the 90% C.L. ellipse of Figure 10.7 in [23], taking the value 1 if on the limiting ellipse. 



T-H± ^ A3 = 0), 
4TO^-t/u^) lines 



White regions are outside the ellipse. The regions inside the black (i) (TO22 = 
magenta (ii) (m^j = ^ A3 = 2m^^^/v'^) and cyan (iii) (m^j = ~"i^± => A3 
fulfill the theoretical constraints for a given 77122 ^^ A3-value. In this plot, A2 = Ai and Ay = Ag. 

orders, p^ will have off-diagonal elements and introduce FCNC, but these will be two-loop 
suppressed [10]. In general one could consider constraints from the renormalization group 
evolution of the Yukawa couplings along the lines of Ref. [15]. 



3 Constraints on the model 

Let us now discuss the constraints that apply to the model. These are analyzed using the 
two-Higgs doublet model calculator 2hdmc [16, 17], where we have implemented the SDM 
as a new model. This makes it straightforward to obtain theoretical constraints, oblique 
parameters, branching ratios^, and cross sections. 

The parameters of the potential are constrained by demanding that the potential be 
bounded from below [7, 18] and that it respects perturbativity and tree-level unitarity[19- 
22]. We will refer to these conditions as "theoretical constraints." We do not list their 
explicit expressions here (see e.g. ref. [16]), but will take them into account in all our 
calculations by using 2hdmc. For details, see ref. [10]. 

The next set of constraints is given by electroweak precision tests (EWPT) , most im- 
portantly those constraints imposed by loop contributions from the scalars to the gauge 
boson vacuum polarizations. Such corrections can be parametrized by the oblique param- 
eters S, T and U [24]. These do not depend explicitly on the potential parameters, but 



^Some of the important decay channels of H and A'^ are one- loop processes that are not included in 
2HDMC, and will be briefly discussed below. They will be calculated and further discussed in [10]. 



do so implicitly through the masses of the scalar particles in the model. All the scalars 
contribute to the oblique parameters, and the resulting expressions are lengthy but not 
particularly illuminating so we do not list them here. We use 2hdmc to evaluate them and 
require the obtained values of S and T to fall within the 90% C.L. ellipse in Figure 10.7 of 
[23]. 

Finally, there are constraints from collider searches at LEP, the Tevatron and the LHC. 
Constraints from LEP and the Tevatron and the 7 TeV LHC constraints (as implemented 
in the latest version of HiggsBounds, v3.8) will in the following be included through 
the use of the HiGGSBoUNDS program [25-27] hnked to 2hdmc. These will be further 
discussed in the next section, but since the A^ and H in our model do not have the same 
interactions as in general 2HDMs, they are not much constrained by previous searches. 
The interactions of the CP-even scalar that has not been observed at the LHC, however, 
include the interactions of the SM Higgs boson, and it could therefore in principle have 
showed up in LHC searches. 

4 Higgs and the LHC 

Let us refer to the particle discovered by ATLAS and CMS as H and to the Higgs boson 
of the SM as .ffsM- In our model, the Ti particle with mass roughly 125 GeV can be either 
the lighter h^, with a heavier H^ remaining to be discovered, or H can be the heavier 
H^, while the h^ is lighter and was not discovered at LEP or the Tevatron due to small 
couplings. We will refer to these two cases as Case 1 and Case 2, respectively. 

As we shall see below, if we want one of our CP-even scalars to account for the recent 
ATLAS data on the ?^ — )• 77 decay channel, we must take the mixing sin a relatively 
large in Case 1 and relatively small in Case 2. In figure 1 we therefore present examples of 
regions in the {mfj± , ?7iA)-plane that are allowed by theoretical and electroweak constraints 
for Case 1 and Case 2. For Case 1 we choose Sa = 0.9, and for Case 2, Sa = 0.1. We 
further use A2 = Ai and A7 = Ag to satisfy eq. (2.12) in the simplest way. Note that 
since the S and T parameters do not explicitly depend on the Aj, only some combinations 
of these parameters (i.e. the couplings and the masses of the scalars) are constrained by 
the electroweak precision data. We therefore additionally show in figure 1 boundaries of 
regions admitted by the theoretical constraints for three different values of A3. 

From figure 1 it is clear that the masses of the scalars should not exceed roughly 
700 GeV in order to fulfill the theoretical constraints. Moreover, the largest allowed region 
is obtained for m\2 ~ 0. In order to give small enough contributions to the S and T 
parameters, mu± and niA must satisfy some approximate custodial symmetries: Define 
M^ = m?^ SIT? a + m^ cos^ a [28]. Then niA ~ 'm^± + 50 GeV when m.^± < M^ or rra^ ~ 
mij± when m?fj± > M^ are allowed, respectively. When m?^^ ^ M^, then < tua < 700 
GeV is allowed. These conclusions are similar to what was found in Ref. [28]. 

In figure 1 we have not yet included collider constraints. Let us estimate the LHC 
constraints on Case 1. Then 7i = h^, while H^ is heavier. It is constrained by the LHC 
searches for the SM Higgs mainly through its decays H^ — > ZZ^*^ — ^ 4^, H^ — )■ l/^VF^*) and 
H^ — > 77. Events where the H^ is produced and decays in one of these channels will look 
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Figure 2. Branching ratios of the i/" boson as a function of its mass. Parameter values are 
TTih — 125 GeV, m/f± = 125 GeV, iua = 175 GeV, A3 — 0, sin a — 0.9, A2 = Ai and A7 — Ag. 



exactly like the corresponding events for a SM Higgs at the same mass, so all experimental 
acceptances are identical. We can therefore compute the a x BR for the relevant channels 
and use the LHC results to constrain the production of H^ . 

The most sensitive exclusion is production via gg — )• H^ and decay via H^ — )• ZZ. 
Let us consider the most recent ATLAS exclusion in this channel shown in Figure 12a of 
ref. [4]. The SM curve in this plot shows a x BR for gg — )■ H^m — ^ ZZ — )■ M. We want 
to rescale this curve to our model. The production processes are the same as in the SM, 
but the Yukawa coupling of H^ to the heavy quark in the loop is suppressed by a factor 
Ca, and the same factor Ca applies to vector-boson fusion and associated production. The 
decay vertex H^ — >• ZZ is also suppressed by the same factor Ca- However, the dominating 
decays to SM particles that contribute to the width T^o are H^ — > ZZ, WW, bb, ti. Again, 
all of these processes are suppressed by the same factor Ca relative to the SM. In addition, 
there may be decays to the new scalars, H'^ -^ h'^h'^ , H+ H~ , A'^A'^ , A'^ Z, H'^W^ , if these 
channels are open. We may thus get a conservative upper bound on the branching ratio 
BR(i:/^^ — )• ZZ) by considering the case where none of the latter channels are open. The 
branching ratio is then the same as in the SM, but if any of the new channels open, it 
becomes smaller. 

The upshot is that to obtain an absolute upper limit on o" x BR in our model, we need 
only rescale the SM result by the factor c^. As an example, in figure la, we have chosen 
Sa = 0.9, so that c^ = 0.19. Referring to the SM curve in Figure 12a of ref. [4], we see 
that rescaling this curve by a factor 0.19, we are below or close to the exclusion limit for 
all masses above 200 GeV. If we instead take Sa = 0.8, our conservative estimate is closer 
to the exclusion region. Note that this is the upper limit in our model, and in most of 
parameter space one or more of the additional decay channels will be open, leading to a 
smaller BR(i?'^ — t- ZZ). We conclude that at least for Sa = 0.9 or larger, a H^ in the mass 
range 200-500 GeV would not have been discovered at the LHC, and this also holds for 
smaller Sa if additional decay channels are open. In this case, the limits become model 
dependent, and we leave this for future studies. As a first illustration of possible additional 
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Figure 3. The maximum fi-Hj-y as described in the text, for (a) Case 1 and (b) Case 2, including 
experimental and theoretical constraints. 



decay channels we calculate the branching ratios for H^ for one point in parameter space, 
as shown in figure 2. It is clear that as soon as the H^ — t- H W^ and/or A^Z channels 
open, they will quickly dominate over the WW and ZZ channels, making H^ more difficult 
to detect. 

We are now in a position to consider the actually observed Higgs boson T-L in our model. 
Similar studies have been done of the Inert Doublet Model in refs. [29, 30]. 

The most significant channel observed is the 77 decay, where the ATLAS experiment 
has observed a slight excess in signal strength, /i = 1.65 it 0.24(stat)_Q']^g(syst), where // is 
the observed signal strength relative to the SM, at a mass of 126.8 GeV [3]. The overall 
best fit mass from all channels is 125.5 GeV. The CMS experiment, on the other hand, sees 
no excess, and has a combined best fit mass of 125.8 GeV [5]. 

The signal strength for 7^ — )• 77 with Ti. = h^ , H^, relative to the SM, is computed as 



^J'H■y^ 



E.<^»(PP^^)BR(^^77) 
Ei (^iiPP -^ Hsm)BR{Hsm -^ 77) 



(4.1) 



where the sums run over all contributing production channels. The production cross sec- 
tions for gg ^ v., vector-boson fusion, and associated production with a vector boson are 
all scaled in the same way compared to the SM: there is a suppression factor s^ for h^ and 
c^ for H^. The signal strength for T-L = h? can then be written 
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(4.2) 



and fJ-HO-yy is obtained in the same way, but replacing s^ ^>- c^. Unlike the case for the 
tree-level decays of the H^ discussed above, the branching ratios here are not equal in 
the SM and in our model, as there are charged scalars running in the loops. We have 
performed a scan over A3 and A7, with —5 < A3 < 5 and —5 < A7 < min{5,A™^^}, and 
in figure 3 plot the maximum obtained signal strength provided that the theoretical, LEP, 
Tevatron and 7 TeV LHC constraints are fulfilled. We plot this in the {mjj±,Sa) plane 



for Cases 1 and 2 discussed above, fixing rriH = 300 GeV for Case 1 and rrih = 75 GeV 
for Case 2. The remaining parameter to be fixed is ttia- In order to fulfill the EWPT 
constraints discussed above, we take, for Case 1, tua = "i_f/± + 50 GeV if wz|^± < M^ or 
niA = fnH± if "^H± ^ ^^' respectively (where M^ was defined above), and for Case 2, we 
fix niA = 'rnjj± . If rufi in Case 2 is increased, the sharp limit at Sq, ~ 0.2 is shifted upwards. 
We have not included off-shell decays Ti -^ //+(*)//-(*) or ^o(*)A°(*) since the widths of 
the H and A^ are very small. 

From figure 3 we see that in order to have a maximal signal strength of one or larger, in 
Case 1, we need Sa ^ 0.8, and in Case 2 we need Sa < 0.2. The enhancements are driven by 
light H^ in the loop, with couplings to the CP-even scalars gtH+H- = —'^^{—SaXs + CaXj) 
and Qhh+h- = — i^lcoAa + SoAy). If the charged scalar in Case 1 is very light (< 80 GeV), 
then it is possible to have almost any value for the mixing down to Sa ~ 0.2 in order to 
obtain the observed 77 signal strength. However, as discussed above, the signal strengths 
for 7i -^ VV, ff are given by ^ihoyv = l^hOff = ^a and Hhovv = f^ROff = c^ if mH± > 
m-^/2. Therefore, from the observations of H ^ ZZ and Ti — >• WW at LHC, the low Sa 
region in figure 3a is disfavored. 

Let us now comment on the phenomenology of the charged scalar. LEP has set limits on 
the mass: the model-independent limit from the Z width is mff± > 39.6 GeV [31], and there 
is a limit of mH± > 80 GeV [31-33] under the assumption that BR(i7+ -^ r+i/)+BR(F+ -^ 
cs) = 1. The charged scalar in our model has no tree- level fermion couplings, and does 
not primarily decay into tu or tb as in the MSSM. The available tree-level couplings are to 
gauge bosons and scalars, but the decays H^ — )• W^^ and H^ — )• W^Z are not allowed at 
the tree-level due to gauge and isospin invariance, respectively. The only possible tree-level 
decays are therefore to final states including at least one other scalar, which can be possibly 
off-shell leading to many-body final states. The usual decay channels tu and tb can indeed 
be generated at the one-loop level, but then we must also consider final states such as VF 7 
and W Z. In [10] we will perform a detailed study of these issues, but we note here that 
only the model independent LEP bound of 39.6 GeV applies, and that the H — t- VF 7 
decay plays an important role, leading to a very interesting phenomenology. 

5 Conclusions 

In the Stealth Doublet Model presented in this paper, the A^ and H have no tree-level 
Yukawa couplings to fermions, so all such couplings are loop-suppressed. Flavor and LEP 
constraints therefore do not apply, and A^ and H^ can be lighter than in standard 2HDMs. 

The CP-even scalars h^ and H^ on the other hand, have very similar couplings to 
standard 2HDMs. We find that either the h^ or the H^ can play the role of the Higgs 
boson that has been observed at the LHC, and refer to these cases as Case 1 for h^ and 
Case 2 for H'^. In both cases there are large allowed regions of parameter space where the 
observed 7^ — )• 77 signal strength is obtained. 

The phenomenology of the model depends on the mixing of the CP-even scalars, given 
by the parameter «„, which also controls the amount of soft breaking of the Z2 symmetry; 
if Sq — 7- or 1, we recover the Inert Doublet Model. In the two cases considered, when 
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the observed 71 = hP , the data imply Sa > 0.8, and when % = H^ , Sa < 0.2 (vahd 
for m/i = 75 GeV). The masses of the charged and CP-odd scalars are constrained from 
the obhque parameters {mfj± ss ruA or m?^± ^ "t-^c^ + i^'jiSa) and tree- level unitarity 
{mfj±,mA < 700 GeV). Moreover, mH± ^ rn-}i/2 is needed to obtain the observed 77 
signal strength. 

It is interesting to note that the parameter space that allows a signal strength for 
7^ — )• 77 of the same magnitude as the one observed at the LHC, is also the same parameter 
space {sa ^ 0.8) where the H^ of Case 1 is safe from having been excluded. If Case 1 of 
our model is realized in Nature, we can therefore expect a discovery of the H^ boson in the 
H^ — )• WW, ZZ, H W^ or A^Z decay channels with more collected luminosity. If Case 2 
is realized, we must instead look for a lighter h^ boson. We will return to these issues and 
the study of the H^ and A^ in a forthcoming publication [10]. 
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